Nearest neighbor exchange in Co- and Mn-doped ZnO 
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We calculate the magnetic interactions between two nearest neighbor substitutional magnetic 
ions (Co or Mn) in ZnO by means of density functional theory and compare it with the available 
experimental data. Using the local spin density approximation we find a coexistence of ferro- and 
antiferromagnetic couplings for ZnO:Co, in contrast to experiment. For ZnO:Mn both couplings 
are antiferromagnetic but deviate quantitatively from measurement. That points to the necessity 
to account better for the strong electron correlation at the transition ion site which we have done 
by applying the LSDA+C/ method. We show that we have to distinguish two different nearest 
neighbor exchange integrals for the two systems in question which are all antiferromagnetic with 
values between -1.0 and -2.0 meV in reasonable agreement with experiment. 
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I. INTRODUCTION 

The manipulation of the electronic spin for infor- 
mation processing gives rise to many advantages and 
would open the way to new applications. This field of 
spintronics has need of a semiconducting, ferromagnetic 
material at room temperature. In that respect, sev- 
eral observations of room temperature ferromagnetism 
in ZnO:Co«iii2, or ZnOiMn'^ had been reported. These 
experimental findings were partly based on theoretical 
predictions using density functional calculations in the 
local spin density approximation (LSDA).'^'* However, 
they are not at all confirmed. Whereas the ferromagnetic 
thin films were mainly produced by laser ablationjSai or 
by the sol-gel method, ^other samples fabricated by pre- 
cursor deposition,'^ molecular beam epitaxy^ or powder 
samples^, showed no signs of ferromagnetism. In that 
contradictory situation we propose to study very care- 
fully the magnetic interaction between two nearest neigh- 
bor substitutional magnetic ions (Co or Mn) in the di- 
lute limit without any codoping effect. Namely, the nu- 
merical values of these exchange couplings are already 
quite well known by susceptibility^ or magnetization step 
measurement^ which allow a careful check of the LSDA 
results. Therefore, we present accurate full potential 
band structure calculations within density functional the- 
ory. We find a rather remarkable discrepancy between 
the measured data and the LSDA results which can be 
considerably reduced by taking into account more prop- 
erly the correlation effects of the transition metal ions 
within the LSDA+f7 method. 

First LSDA studies found ZnO: Co to be ferromag- 
netic, but ZnO:Mn antiferromagnetic* Pseudopotential 
calculations on large supercclls^'^^ which were performed 
later on, could specify the different couplings more in de- 
tail. They found a competition between ferromagnetic 
and antiferromagnetic interactions in ZnO:CG^^*' which 
we confirm and argued for the necessity of additional 
electror&iS or hole doping^"' to stabilize ferromagnetic 



order. (Such an additional doping will not be studied 
here.) But the presence of ferromagnetic and antifer- 
romagnetic couplings at the same time would lead to a 
very small Curie- Weiss constant in contrast to the ob- 
served one which is clearly antiferromagnetiCifi'S- - Even 
more clear is the contradiction between the measured 
data and the LSDA results for ZnO:Mn: magnetiza- 
tion step measurements^ lead to antiferromagnetic ex- 
change integrals of -1.56 and -2.08 meV for the two near- 
est neighbor positions possible, although the experiment 
does not allow to assign these values to certain bonds 
in an unambiguous way. Transforming the energy differ- 
ences into numerical values of exchange integrals (which 
was not done in Refs. 9 10) we will show that LSDA over- 
estimates the experimental results considerably. After 
having stated the discrepancy between LSDA and ex- 
periment, we show that the LSDA-|~J7 method can cure 
these deficiencies and leads to reasonable exchange cou- 
plings. The importance of the LSDA-|-t/ approximation 
has also recently been pointed out for the related system 
ZnSe:Mnii 



II. SUPERCELL CALCULATIONS 

To determine the nearest neighbor exchange couplings 
we performed several supercell calculations. ZnO crys- 
tallizes in the hexagonal wurtzite structure (space group 
P63mc) with the lattice parameters a — 3.2427 A and 
c = 5.1948 A.^'^ We consider here pure substitutional de- 
fects and neglect the influence of lattice relaxations. Due 
to the wurtzite structure there are two crystallograph- 
ically different nearest neighbor positions: the in-plane 
nearest neighbor within the plane perpendicular to the 
hexagonal axis c and the out-of-plane nearest neighbor. 
Their magnetic couplings were studied by using 4 differ- 
ent supercells (A,B,C, and D). Each supercell is formed 
by multiples of the primitive lattice vectors a, b, and 
c, like the 2x2x1 supercell A, shown in Fig. 1. The 
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FIG. 1: (Color online) Crystal structure of the supercell A. 
(TM = Co, Mn) 
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FIG. 2: (Color online) (a,6) plane corresponding to supercells 
A and C as well as (|(a + 6), c) plane corresponding to super- 
cells B and D, with the definition of in-plane and out-of-plane 
exchange constants Jin and Jout, respectively. (TM — Co, 
Mn) 



supercells A and C probe the in-plane nearest neighbor 
exchange Jin (see Fig. by a chain of Co (or Mn) im- 
purities (supercell A) or an isolated pair (3x2x1 super- 
cell C). The out-of-plane exchange Jout is probed by the 
2x2x1 supercell B (chain) and the 2x2x2 supercell D 
(pair) . The so defined supercells A-D coincide with those 
used in Ref. jS- 

The supercell calculations were performed using 
the full-potential local-orbital (FPLO) band structure 
scheme,^^ In the FPLO method a minimum basis ap- 
proach with optimized local orbitals is employed, which 
allows for accurate and efficient total energy calculations. 
For the present calculations we used the following basis 
set: Zn, Co(Mn) 3s3p-As4p3d, O 2s2p;3d. The inclusion 
of the transition metal 3s and 3p semicore states into 



the valence was necessary to account for non-negligible 
core-core overlap, and the O 3d states were used to im- 
prove the completeness of the basis set. The site-centered 
potentials and densities were expanded in spherical har- 
monic contributions up to Imax = 12. 

The exchange and correlation potential was treated 
in two different ways. First, the local spin-density ap- 
proximation (LSDA) was used in the parametrization of 
Perdew and Wangii^ However, as will be shown more in 
detail below, this approximation has severe deficiencies 
in the present case. The Co(Mn) 3d states are in real- 
ity more localized than in the LSDA calculation. This 
correlation effect was taken into account by using the 
FPLO implementation of the LSDA-t-[/ method in the 
atomic limit schemeiiSiii The convergence of the total 
energies with respect to the k-space integrations were 
checked for each of the supercells independently. We 
found that 8x8x8 = 512 k-points were sufficient in 
all cases, and this parameter was used in the calculations 
reported below. 



III. EXCHANGE COUPLINGS 

The Zn^+ ion in ZnO has a completely filled 3d shell 
and, correspondingly, no magnetic moment. If Zn is 
replaced by Co or Mn, the valence 2+ is not changed, 
which is also proved by our bandstructure results below. 
It means that these substitutional impurities provide no 
charge carriers. The configuration of Co^+ is and that 
of Mn^+ is d^. Therefore, they have a spin S = 3/2 or 
S = 5/2, correspondingly. The tetrahedral and trigonal 
crystal fields, together with the spin-orbit coupling, lead 
to a magnetic anisotropy. But they are not large enough 
to destabilize the high-spin states in the given cases, as 
supported by electron paramagnetic resonance and mag- 
netization measurements (see Refs. I^18lll9l for ZnO:Co 
and Ref. ia for ZnO:Mn). 

The present work is devoted to determine the domi- 
nant exchange couplings between two localized magnetic 
ions. It can be expected that the dominant couplings oc- 
cur between nearest neighbor impurities, each carrying a 
local spin S^. Then, the Heisenberg Hamiltonian for a 
localized pair of spins is given by 



H — —2JSiSj 



(1) 



The corresponding total energies for ferromagnetic (FM) 
and antiferromagnetic (AFM) arrangements of the two 
spins are 

Efm = -J [St{St + 1) - 25(5 + 1)] , 
Eafm = J [25(5 + 1)] , (2) 

with the total spin St — 25 of two parallel spins 5. This 
leads to the energy difference between the FM and AFM 
states per magnetic ion: 
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with = 3 or 5 for Co or Mn. That energy difFer- 
ence can be compared with the corresponding energy 
differences of isolated pairs in the supercells C and D. 
The supercells A and B, however, correspond to chains 
of magnetic ions. For our purpose, it is sufficient to use 
an approximative expression for the energy difference be- 
tween FM and AFM states of a Heisenberg chain. Then, 
each magnetic ion has two nearest neighbor magnetic ions 
which doubles the previous energy difference P|l: 
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That result can also be calculated more explicitly by 
comparing the energies of the FM and AFM states of 
a Heisenberg chain decomposed into a series of pairs. 

It is remarkable that the so defined exchange cou- 
plings are experimentally measurable. They will be de- 
noted here as J™ and Jout for the in-plane and out-of- 
plane nearest neighbors, respectively. The most precise 
measurements had been performed using magnetization 
stepa^ for the case of ZnO:Mn which leads to different 
values for Jin and Jout ■ We are not aware of such mea- 
surements for ZnO:Co. An average value of Ji„ and Jout 
is however accessible by susceptibility measurements i 



IV. RESULTS 
A. ZnO:Co 

The density of states (DOS) of the FM solution (for 
supercells A and B) is shown in Fig.O Its main features 
agree with previous calculations >Si The minority and ma- 
jority Co 3d states are found to lie mainly in the gap 
between the valence band of predominantly oxygen 2p 
character and the Zn 4s-4p conduction band. The Zn 3d 
states are located at about -7 eV at the bottom of the va- 
lence band, deep below the Fermi level. The occupation 
of the Co 3d level is close to Sd^ and the magnetiza- 
tion Mf^(Co) = 2.6^B, i.e. rather close to S* = 3/2. 
The majority Co 3d states are located just above the 
oxygen valence band. There is a small hybridization of 
the minority 3d level with the conduction band which 
makes the material half-metallic in the LSDA (see also 
Refs.'^'T^. However, this half- metallic character would 
correspond to a partial electron doping and is an artifact 
of the LSDA solution. 

Already the gap of pure ZnO (experimental gap: 3.3 
eV) is underestimated by LSDA (FFLO leads to a gap of 
1.4 eV). To study the 3d levels of an isolated Co impu- 
rity we calculated the electronic structure of a CoZn708 
supercell by LSDA. The corresponding DOS of the FM 
solution is rather close to that one shown in Fig. O but 
the band structure (not shown) allows a better analysis 
of the impurity levels, free of hybridization effects be- 
tween neighboring impurities. For instance, the crystal 
field splitting of the Co 3d levels is clearly visible at the 
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FIG. 3: (Color online) Comparison of total and partial DOS 
(for 1 Co atom) for supercells A and B, ZnO:Co, in the FM 
case obtained by LSDA. 



F-point. In each of the spin channels the two- fold degen- 
erated Cg levels are situated below the t2g orbitals, for 
majority spin by about 0.5 eV. The eg levels remain de- 
generate in the trigonal case, whereas the t2g ones split 
into a lower singlet and an upper doublet. This trigonal 
splitting is of the order of 0.2 eV. Due to hybridization 
effects, the Co 3c? impurity band is considerably broader 
in Fig. 121 than for CoZnrOs by roughly a factor of three 
for majority spin. That shows that hybridization and 
crystal field splitting are of the same order of magnitude 
for two neighboring Co impurities in ZnO. 

The differences of total energies between FM and 
AFM solutions give the corresponding exchange cou- 
plings in the way described above. The values for in- 
plane and out-of-plane exchange Jm and Jout are col- 
lected in Tables I and II. Like in the previous pseudopo- 
tential calculation32ii£ the in-plane exchange is antiferro- 
magnetic, but the out-of-plane exchange ferromagnetic. 
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TABLE I: Calculated in-plane exchange Jin for ZnO:Co using 
LSDA {U = 0) and LSDA-Kt/ {F° = U ^ 0, = 7.9 eV and 
— 5.0 eV) for the supercells A and C. Also given are the 
corresponding energy differences per Co ion and the magnetic 
moments. 
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TABLE II: Calculated out-of-plane exchange Jout for ZnO:Co 
using LSDA {U = 0) and LSDA+f/ {F° ^ U ^ 0, F'^ = 
7.9 eV and = 5.0 eV) for the superceUs B and D. Also 
given are the corresponding energy differences per Co ion and 
the magnetic moments. 



The differences between supercells A and C (or between 
B and D, correspondingly) arise due to finite size effects 
or deviations from the Heisenberg model. Qualitatively, 
all available LSDA energy differences agree among each 
other and with our FPLO results (see Table III). But 
there are rather remarkable numerical deviations between 
the different methods. Please note, that the previous 
authors^- did not convert the LSDA energy differences 
into exchange couplings. That we have done to allow the 
comparison with experimental data. 

The competition between FM and AFM nearest neigh- 
bor exchange in ZnO:Co is in contrast to experimen- 
tal results^'^'^ which show dominantly AFM couplings. 
Other problems of the LSDA solution are the following: 
(i) the semi-metallic character, (ii) the insufficient local- 
ization of the Co 3d states, and (iii) the position of the 
impurity 3c? levels. Namely, photoemission spectroscopy 
shows them as deep impurity levels close to the valence 
bandit But the experimental energy difference to the 
top of the valence band of only 0.4 eV is smaller than the 
corresponding distance of the center of gravity of the Co 
3d level (about 1 eV for majority spin). 

So, we should look for a theoretical method which takes 
into account correlation effects more properly. One has 
to distinguish the correlation in valence and conduction 
band leading to the incorrect gap value and the correla- 
tion effects in the Co 3d orbitals. The first effect might 
be repaired by the GW approximation,^^ or including 
the self interaction correction (SIC) as proposed in Ref. 

But we do not expect that it would considerably 
improve the exchange couplings. The correlation in the 
Co 3d shell will be taken into account in our present 
work by the LSDA-f[/ scheme using the atomic limit 
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FIG. 4: (Color online) DOS of ZnO:Co, superceU A, in the 
FM case, calculated with LSDA-I-?/. 



functionaliSi We also tried the "around mean field" ver- 
sion of LSDA-f [/ which gives similar results as reported 
below for ZnO:Co, but which leads to no improvements 
for ZnO:Mn due to the known peculiarities of this ftmc- 
tional for the d^ configuration. ^^'-^ The parameters for 
ZnO:Co have been chosen similar to those for CoO^ 
namely the Slater parameters and such that the 
Hund's rule exchange j" = (i^^ _^i^'')/14 = 0.92 eV and 
the ratio F^/F^ = 0.625 are close to ionic values2^ which 
leads to = 7.9 eV and = 5.0 eV. The parameter 
F'^ = U is less well known since it is more affected by 
screening effects. We found that a value of at least 3 eV 
is necessary to stabilize an insulating solution. In the re- 
gion of realistic F*^ parameters its influence on the values 
Jin and Jo„t is small as indicated by the comparison of 
results for F° = 6 and 8 eV in Tables I and II. 

The DOS (Fig. ^ shows then clearly an insulating 
state and the occupied Co levels are much closer to the 
valence band than in LSDA in better agreement with 
photoemission dataiSi The gap value {Eg « 0.7 eV for 
the FM case, F° = 6 eV) is not much improved, but 
that can also not be expected since we did not change 
the potentials for oxygen or zinc s-p states. However, 
it is remarkable that now both exchange couplings are 
antiferromagnetic. 

We find the in-plane exchange Jm = — 2.0 meV to be 
larger than the out-of-plane exchange Jout = — 1.0 meV 
(for U — 6 eV). Experimentally, the antiferromagnetic 
nearest neighbor exchange was determined to be J = —33 
K or -2.8 meV from the high-temperature Curie- Weiss 
constant of the magnetic susceptibility^ which exceeds 
slightly our values. 
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TABLE III: Available theoretical LSDA results for Jin and 
Jout (all data in meV) of ZnO:Co and ZnO:Mn taken from the 
literature. The published energy differences were converted 
into exchange constants. 



B. ZnOrMn 

In Fig. 5, we show the DOS of ZnO:Mn, calculated 
with supercell A and a ferromagnetic arrangement of the 
Mn moments. As in the case of ZnO:Co, LSDA yields 
a metallic solution. The Mn Zd shell is approximately 
half filled and the exchange splitting between the centers 
of gravity of the occupied majority and unoccupied mi- 
nority subbands is about 3.5 eV. The total spin moment 
4.96/XB/Mn is close to the expected S = 5/2 value. The 
Mn atoms carry a spin moment of about but also 

their four nearest neighbor O atoms have a weak induced 
spin moment. The Mn 3d impurity states are mainly lo- 
cated in the upper part of the ZnO gap and weakly hy- 
bridize with the Zn 4s-4p conduction band, which gives 
the solution a metallic character. Measurements of the 
band gap of ZnO:Mn films^^ (see Ref. 27 for an overview) 
find a slight blue shift of the absorption edge with a signif- 
icant amount of mid gap absorption above 2.5 eV. This is 
consistent with a position of the impurity levels around 
the upper edge of the valence band in contrast to the 
LSDA result. 

As shown in the lower graph of Fig. 5, LSDA+t/ shifts 
the highest occupied Mn 3d levels to the top of valence 
band, so that the solution becomes insulating. The pa- 
rameters used in the calculation are U ^ F'^ ~ 6 eV, 
= 7.4 eV and = 4.6 eV, corresponding to = 
0.86 eV, the value chosen for MnO in Refs. 17 25. As in 
the case of ZnO:Co, the value of the band gap Eg « 0.4 
eV is smaller than the experimental one (see discussion 
above), but the position of the Mn 3d impurity levels is 
considerably improved. Compared to the LSDA calcula- 
tion, the partial Mn 3d DOS is slightly broadened and the 
unoccupied Mn 3d minority spin states are shifted fur- 
ther away from the Fermi level. The total spin moment 
is now 5/iB corresponding to an ideal S — 5/2 situation 
with the magnetic contributions almost entirely due to 
Mn (see Tables IV and V). 

In contrast to the case of ZnO:Co, LSDA yields an 
AFM exchange coupling for both types of nearest neigh- 
bor pairs in ZnO:Mn. This is in qualitative agreement 
with the magnetization step measurements of Ref. IH 
where values Jin = —2.08 meV and Jout = —1.56 
meV have been obtained. However, the LSDA values 
Jin = —4.9 meV and Jout = —4.1 meV are 2-3 times 
larger than the experimental ones. Similar results have 
also been obtained by Sluiter et aliM, who find both cou- 
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FIG. 5: (Color online) DOS of ZnO:Mn, comparison of LSDA 
and LSDA+f/, supercell A. 



plings to be strongly AFM (see Tab. III). However, the 
rather poor quantitative agreement between the calcu- 
lated and measured J values indicates that the Mn im- 
purity levels are not well described within LSDA. Again, 
the LSDA-I-C/ functional strongly improves the agree- 
ment of the calculations with experimental data. Taking 
U = 6 eV, we obtain Jin — —2.0 meV and Jout — —1.3 
meV, close to the measured values. As can be seen in Ta- 
bles IV and V, the calculated exchange couplings depend 
only moderately on the choice of the J7-parameter. The 
larger coupling is obtained for the in-plane pairs despite 
the larger distance, as was already assumed in Ref. [ 



V. DISCUSSION 

The reason for the competition of FM and AFM ex- 
change couplings within LSDA for ZnO: Co is schemati- 
cally shown in Fig. Let us first consider an isolated 



6 



Cell 


U[eV] 


A E[meV/Mn] 


Jin [mcV] 




Mf*^ (Mn) 


A 





147 


-4.9 


4.52 


4.62 


C 





74 


-4.9 


4.57 


4.62 


A 


6 


59 


-2.0 


4.87 


4.89 


A 


8 


48 


-1.6 


4.94 


4.94 



.4ITVI Couiilhig; 



FM C'ouijliiig 



TABLE IV: Calculated in- plane exchange Jin for ZnO:Mn 
using LSDA {U = 0) and LSDA+(/ (F° = U ^ 0, f ^ = 
7.4 eV and = 4.6 eV) for the supercells A and C. Also 
given are the corresponding energy differences per Mn ion 
and the magnetic moments. 
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TABLE V: Calculated out-of plane exchange Jout for ZnO:Mn 
using LSDA {U = 0) and LSDA+(7 {F° ^ U ^ 0, = 
7.4 eV and F"^ = 4.6 eV) for the supercells B and D. Also 
given are the corresponding energy differences per Mn ion 
and the magnetic moments. 



Co ion in ZnO. The 3c? levels are split by the crystal 
field (CF) into lower Cg and upper t2g levels. They are 
also influenced by the exchange splitting between spin up 
and spin down electrons caused by the local Hund's rule 
exchange. These local energy levels are filled with 7 elec- 
trons in the case of Co. The LSDA-DOS (Figs. 3,5) shows 
that the CF splitting is smaller than the exchange split- 
ting. Fig.lHlpresents the hybridization effect on the Co 3d 
energy levels if the two Co ions come close together. The 
hybridization leads to the formation of a pair of bond- 
ing and antibonding hybrid orbitals for each id energy 
level Ei of Co. The bonding and antibonding orbitals 
have energies Ei — ASj and Ei + AiJ^, correspondingly. 
Therefore, the complete filling of these two orbitals does 
not lead to an energy gain, but a partial filling does. In 
such a way, the energy gain for an AFM arrangement of 
spins is evident. For a FM arrangement, the energy gain 
is only possible by a crossing of the eg and t2g levels for 
minority spin. This competition between the FM and the 
AFM energy gain is apparently not identical for in-plane 
and out-of-plane exchange, leading to different signs of 
the exchange couplings. However, as already discussed, 
that is an artifact of the LSDA solution. 

In LSDA-)-?7 (Fig. [71) the unoccupied minority spin 
energy levels are much higher than the occupied ones. 
Therefore, the crossing of minority Cg and t2g energy lev- 
els, and also the FM energy gain, is not possible. As a 
consequence, one finds an AFM superexchange coupling 
in ZnO:Co independent of the geometrical configuration. 

The situation for ZnO:Mn is different. In that case, 
only the majority spin is completely filled with 5 elec- 
trons and the minority spin is nearly empty in LSDA, 
and completely empty in LSDA-I-C/. Since the exchange 
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FIG. 6: (Color online) Schematic view of hybridization and 
CF effects on two close Co 3ci shells in the LSDA case, com- 
parison between AFM and FM couplings. The energy levels 
of majority spin are lower in energy than minority spin ones 
due to the Hund's coupling J^ . 
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FIG. 7: (Color online) Schematic view of hybridization and 
CF effects on two close Co 3d shells in the LSDA-|-(7 case, 
comparison between AFM and FM couplings. The energy 
levels of majority spin are lower in energy than minority spin 
ones due to the Hund's coupling J^ . The occupied energy 
levels are lower in energy than unoccupied ones due to the 3d 
shell correlation effect U . 



splitting is larger than the CF splitting, there is no en- 
ergy gain possible for a FM arrangement of impurity 
spins neither in LSDA nor in LSDA-f J7. The distance 
between occupied and unoccupied energy levels increases 
in LSDA-I-J7. Therefore, the energy gain is reduced and 
leads to exchange couplings which are much closer to the 
experimental values obtained by magnetization step mea- 
surements than those obtained by LSDA. It should be 
noted that our calculation confirms also the assignment 
of Ref. H that Ji„ corresponds to the largest coupling. 

Of course, there are still numerical error sources on the 
exchange couplings Jin and Jout which were calculated by 
LSDA+C/. First of all, we should note the poor knowl- 
edge of correlation parameters U , F^, and F'^, which 
influences the results. Second, there might still be finite 
size effects due to the specific form of the supercells cho- 
sen. And finally, also a small basis set dependence of 
the FFLO method cannot be excluded. All together, we 
would estimate an upper error of about ±30 per cent for 
the calculated exchange couplings. 
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VI. CONCLUSIONS 

In the dilute limit, nearest neighbor pairs of Co and 
Mn-impurities in ZnO have antiferromagnetic exchange 
couplings. That is the result of theoretical calculations 
which take into account the electron correlations in the 
impurity 3d shell properly, and is in agreement with the 
experimental results. This AFM nearest neighbor ex- 
change excludes ferromagnetism for pure substitutional 
Co or Mn defects in ZnO in the dilute limit. The observed 
FM in ZnO: Co and ZnO:Mn should have a different ori- 
gin. There are several proposals in the literature like sec- 
ondary phases:^ - or cation vacancies or other defectsi22i2£ 

The LSDA predictions might be misleading and should 
be considered with care since they do not correctly take 
into account the localized character of the transition 



metal impurities. On the contrary, the LSDA-I-C/ val- 
ues are in good agreement with experimental exchange 
constants derived from magnetization step measurements 
and high-temperature susceptibility data.^'^ So our study 
puts considerable doubts on the value of pure LSDA pre- 
dictions (as published for instance in Refs. I MTol) . at least 
in the case without additional electron or hole doping. 
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